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A simulation of quenched QCD with the overlap Dirac operator has been completed using 100 Wilson gauge 
configurations at [3 = 6 on an 18 3 x 64 lattice. We present results for meson and baryon masses, meson final state 
"wave functions" and other observables. 



1. Simulation Details 

We have studied quenched QCD on an 18 3 x 64 
lattice, with the Wilson gauge action at f3 = 6 
and with Neuberger's overlap Dirac operator for 
lattice fermions. We used a sample of 100 gauge 
configurations separated by 10,000 upgrades of a 
6-hit Metropolis algorithm (acceptance ~ 0.5). 
A previous preliminary analysis £t] used the first 
60 quark propagators that were produced; results 
presented here include all 100 quark propagators. 
We refer to PP and the references therein for more 
background information and details of the over- 
lap operator used. New results include the vector 
meson spectrum, the use of extended sink opera- 
tors in constructing meson correlation functions, 
and observation of quenched chiral logs. Results 
for kaon weak matrix elements are presented in a 
separate contribution 0. 

From Wilson loops, we obtained r$/a = 5.36 ± 
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0.11 (a- 1 = 2.11 ± 0.04 GeV) for the Sommer 
scale defined by r^Flr^) = 1.65, ro = 0.5 fm. 

We calculated overlap quark propagators 
for a single point source, for all 12 color- 
spin combinations and quark masses am q = 
0.03,0.04,0.06,0.08,0.1,0.25,0.5,0.75, using a 
multi-mass solver. 

The Zolotarcv approximation with 12 poles was 
used for the first 55 configurations. For the re- 
maining 45 configurations, the Chcbyshcv ap- 
proximation (degree 100 ~ 500) was used, after 
Ritz projection of the lowest 12 eigenvectors of 
H 2 , as it was found to provide about 20% better 
performance. The convergence criterion used was 
|Z)t Dip - x\ 2 < 10" 7 . 

The computations were performed with shared 
memory F90 code, optimized and run on 16 and 
32 processor IBM-p690 nodes at BU and NCSA. 



2. Meson Results 

We begin with a verification of the axial Ward 
identity (AWI) prediction for the ratio 



p(t) = G 
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(t)/G PP (t) 
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as presented in Fig.^ The fit to ap = A+a(m ql + 
m q2 )/Z A + C[a(m ql + m q2 )} 2 gives A = (3.09 ± 
9.57) x 10~ 6 , Z A = 1.5570 ± 0.0004, and C = 
(-5.12 ± 0.61) x 10~ 3 , as shown in Fig. H 
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Figure 1. AWI ratio as a function of time for all 
denerate quark mass combinations. 




Figure 3. Time dependence of a pseudoscalar cor- 
relator. Its value on each gauge configuration as 
well as its average value are shown. 
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Figure 2. Axial Ward identity, continued. 
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Figure 4. Scanning the window for optimum t Tl 



2.1. Meson Spectrum 

Pseudoscalar density 2-point functions were 
formed from quark propagators and then fit to 
the usual cosh form to calculate the ground state 
mass. An example correlator is shown in Fig. [3] 
for input quark masses am q i = 0.08 and am q 2 = 
0.10. 

Effective mass plots are often used to deter- 
mine the fitting window {t rn i n ,t ma x} for the cor- 
relation function. In some cases (e.g. for vector 
mesons), determining the best window can be dif- 
ficult. Therefore, we used a method that scanned 
the possible windows, and then chose the smallest 
value of t m i n (consistent with the errors) allowed 
before the clear effect of higher states caused the 
mass prediction to rise. Use of a cosh fit allowed 



tmax to be fixed at 32. An example of such a scan 
is shown in Fig. The pseudoscalar spectrum 
for all possible input quark mass combinations is 
shown in Fig. |S| Using a larger lattice than [3J 
allowed for observation of vector meson states, as 
shown in Fig. H3 

2.2. Extended sink operators 

Since mesons are not point particles, an in- 
creased signal can be obtained by using icx- 
tended (or ismeared) operators at the source and 
sink [3]. Matrix element calculations required 
point sources and generating an additional set 
of propagators with extended sources would have 
been computationally too costly. Thus we limited 
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Figure 5. Pscudoscalar spectrum. 
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Figure 6. Vector meson spectrum for point and 
extended sinks. 
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Figure 7. Extended sink functions C(r, t) for var- 
ious t/a. 
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Figure 8. Vector meson spectrum for am q < 0.1. 



our analysis to extended sinks. 

The use of extended sinks allowed for a larger 
fitting window and improved the prediction for 
mesonic observables. The extended sink correla- 
tion function is 

Gf(t) = (E V^)rv>(2/)/(x - y)VXO)r^(O)) (2) 

where x = (x, t) and y = (y, t), and the gauge 
is fixed to Landau gauge. With /(x — y) = 
S r (\x — y|— r), we calculate the quantity C(r, t) = 
G Sr {t) shown in Fig. We the n use the 
asymptotic value Co(r) = C(r, oo) to define 
an extended sink correlation function G ex t(t) = 
Gc (t) ■ The functions C (r, t) stabilizes at ap- 
proximately t/a = 10, representing the meson fi- 
nal state "wave function." 



The use of extended sinks was most valuable in 
the calculation of the vector meson spectrum, due 
to the increase in the size of the fitting window. 
The point sink and extended sink vector meson 
spectra are compared in Fig.|H|for low quark mass, 
am, < 0.1. A fit of the quark-mass dependence 
of the vector meson mass obtained from the ex- 
tended sink data gives 

aM = 0.4159(11) + 1.044(16)a(m 9 i + m, 2 )/2 .(3) 

2.3. Quenched chiral logs 

The sharp rise at low am q of the ratio of the 
pscudoscalar mass squared (aM) 2 to (am q \ + 
am q 2 ) , which we observe in Fig. [51 can be inter- 
preted as an indication of quenched chiral logs. A 
value for the quenched chiral log parameter 6 was 
obtained via a fit to the expression y — 1 + S x 5 , 
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Figure 9. Evidence for quenched chiral logs. 



Figure 10. Caution about fluctuations. 



where 



2m! Mf 2 
mi + m 2 Mi! 



2m 2 
mi + m 2 



Mf 2 



and 



x = 2+ mi+m 2 frn, 
mi — m 2 V mi 



(4) 



(5) 



The resulting value for 8 in the PP-SS channel is 
r5 = 0.13 ±0.09. 

We fit the degenerate quark mass PP-SS chan- 
nel results to the expression 



(Mo) 2 = A{ma) 1/{1+5) + B{maf 



(6) 



where M is the PP-SS mass and m = m q \ + m q2 . 
The fit is shown in Fig. El an d the results for the 
parameters are 

A = 0.441(1), 5 = 0.604(3), 5 = 0.243(1), (7) 

consistent with the value of 6 obtained above as 
well as elsewhere |5I6I7) . Note that a fit to (Ma) 2 
in Eq. JfjJ) containing an additional constant term 
C does not change the central values of A, B, or 
S and produces a value for C that is very small 
and consistent with zero. 

3. Baryon Results 

In calculating the baryon spectrum, we found 
that caution is needed in determining an appro- 
priate fitting window for the correlation func- 
tions, due to fluctuations at large t. 
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Figure 11. Baryons: octet and decuplet spectra 
for various fitting windows. 



Shown in Fig. ar e the average baryon prop- 
agators for the octet and decuplet, after multi- 
plication by cxp(0.65£/a) and the mapping y — > 
2dtam%/(2d)], d = 10~ 6 . The mapping allows 
for a clear exposition of the fluctuations. 

Fig. ^2 shows the octet and decuplet baryon 
masses vs. total quark mass for all quark mass 
combinations where two of the quarks have the 
same mass m q i, while the third quark has mass 
m,j 2 , and for various correlator fitting windows. 
The sensitivity to the fitting window is appar- 
ent. The ultimate fitting window chosen was 
8 < t/a < 15. 

Fig. 1121 shows the low quark mass spectrum, 
with fitting window 8 < t/a < 15, for the two 
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Figure 12. Spectra for the two octet baryon states 
discussed in the text. 



Figure 13. Baryons for am qi — am q2 < 0.1. 



octet states 

81 = (| TIT ) - 1 ITT ))/V2 (8) 

8 2 = (| TIT) + I ITT)-2| TU))/Ve (9) 

with am q i = 0.03 and 0.03 < am q i < 0.25. The 
corresponding correlation functions become iden- 
tical for m q 2 = m q i, but would describe A, £-like 
states, respectively, for m q \ ^ m q 2. 

Fig. El shows the degenerate quark mass 
octet and decuplet spectrum for low quark mass 
(am q < 0.1) again with fitting window 8 < 
t/a < 15. Extrapolating to the chiral limit we 
get aM 8 = 0.562(4), aM w = 0.747(3), aM p = 
0.4159(11), and correspondingly 

M p /M 8 = 0.740(7) (10) 

M w /M s = 1.329(11) (11) 

while experimentally M p /Mn = 0.820 and 
M A /M N = 1.312. 

4. Conclusion 

Our results show that use of the overlap formu- 
lation is practical, albeit computationally costly, 
for the calculation of quark propagators on large 
lattices, and validate the good chiral properties 
of this discretization. Work in progress includes 
the evaluation of quark correlation functions in 
baryon final states and quenched calculations on 
a 14 3 x 48 lattice at f3 = 5.85 to investigate scal- 
ing. 
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